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• • . 1 INTRODUCTION 



ABSTRACT 

The weak lensing surveys have the potential to probe directly the clustering statistics 
of dark matter in the universe. Recent studies have shown that it is possible to predict 
analytically the whole probability distribution function (pdf) and the bias associated 
with the collapsed objects in the highly non-linear regime using the hierarchical ansatz. 
We extend such studies to the quasi-linear regime where the hierarchical ansatz is re- 
placed by the tree-level perturbative calculations to an arbitrary order. It is shown 
how the generating function techniques can be coupled with the perturbative calcu- 
lations to compute the complete pdf and the bias in the quasi-linear regime for the 
weak-lensing convergence field. We study how these quantities depend on the smooth- 
ing angle and the source red-shift in different realistic cosmological scenarios. We show 
that it is possible to define a reduced convergence whose statistics is similar to un- 
derlying 3D mass distribution for small smoothing angle but it resembles projected 
mass distribution for large smoothing angles. We have also compared our peturbative 
results with log-normal model for pdf and bias and found a good agreement between 
the two analytical results. 
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The theoretical study of the weak gravitational lensing to probe the statistics of large scale structures in the universe has 
received much attention (Mellier 1999, Bartelmann & Schneider 1999) in recent past. 

Based on the earlier works done by several authors (Gunn (1967), Blandford et al. (1991), Miralda-Escude (1991) and 
Kaiser (1992)), current development has focussed mainly on two different directions. Much progress have been achieved in the 
field of numerical simulations of the weak lensing, in particular the construction of shear maps which provide invaluable tools 
in testing the analytical results. The numerical simulations typically employ N-body simulations, through which ray tracing 
experiments are conducted (Schneider & Weiss 1988; Jarosszn'ski et al. 1990; Lee & Paczyn'ski 1990; Jarosszn'ski 1991; Babul 
& Lee 1991; Bartelmann & Schneider 1991, Blandford et al. 1991). Building on the earlier works of Wambsganns et al. (1995, 
1997, 1998) the most detailed numerical study of lensing was done by Wambsganns, Cen & Ostriker (1998). Other recent 
studies using ray tracing experiments have been conducted by Premadi, Martel & Matzner (1998), van Waerbeke, Bernardeau 
& Mellier (1998), Bartelmann et al (1998), Couchman, Barber & Thomas (1998) and Jain, Seljak & White 2000. 

On the analytical front, Villumsen (1996), Stebbins (1996), Bernardeau et al. (1997) and Kaiser (1992) have focussed on 
the use of perturbative techniques to study weak lensing. More recently it has been shown that the hierarchical ansatz for the 
higher order correlation functions can be used to describe the statistical properties of the weak lensing convergence fields for 
the case of small smoothing angles (Hui 1999, Munshi & Coles 1999, Munshi & Jain 1999, Valageas 1999a,b, Munshi 2000, 
Munshi & Coles 2000a,b). 

In this paper we will show that while perturbative calculations have already been used to compute the lower order 
moments of convergence field, (see e.g. Bernardeau, van Waerbeke & Mellier 1997, van Waerbeke, Bernardeau & Mellier 1998, 
for detailed review on perturbative calculations in the context of weak lensing surveys see Bernardeau 1999) ; when coupled 
with generating function techniques it can also reproduce the complete pdf and bias functions associated with the smoothed 
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projected density field. We also show that these results can simply be obtained by replacing the generating function of the tree 
hierarchy in non- linear regime with its quasi- linear counterpart. While the non- linear generating function corresponds to the 
three dimensional generating function for the matter correlation hierarchy, in the quasi-linear regime it has to be replaced by 
a generating function which represent the smoothed two-dimensional projected density field. These generating functions have 
already been studied in great detail by Bernardeau (1992, 1994, 1995, 1996) for the case of top-hat smoothing function. Using 
the recent analytical results from the study of statistics of the convergence field for small smoothing angles, we will introduce 
a reduced convergence field r/ and study how its statistical properties (including the lower order moments) are related with 
its counterpart used to describe the projected density field. 



2 THE STATISTICS OF THE CONVERGENCE FIELD FOR LARGE SMOOTHING ANGLES 

The weak lensing convergence k is simply the projected mass distribution and its statistics carries valuable information about 
the geometry and dynamics of the background universe as well as the physics of collision-less clustering. 

k( 7 ) = / d X Lo{ X )5{r{x)~l) (1) 
Jo 

Throughout our discussion we will be placing the sources at a fixed red-shift (an approximation not too difficult to modify 
for more realistic description). The weight function can be expressed as lo(x) = 3/2ac~~ 2 Hq Q, m r(x) r (Xs — X)/ r {Xs)- Where 
Xs is the comoving radial distance to the source placed at a red-shift z s . We will be focusing on the quasi-linear regime where 
it is possible to expand the density contrast using a perturbative series. 

K M (7 ) + K P) (7) + . . . = d X oj{x){5 (1) {r{xh) + 5 (2) (r( X )7) + • • •) (2) 
Jo 

It will be useful to work in Fourier domain and the Fourier decomposition of S can be written as: 

«(7)= / dxco{x) / — —3 exp(ixfc|| + ir6k±)5 k (3) 
Jo J ( 27r ) 

Where we have used k\\ and k± to denote the components of wave vector k, parallel and perpendicular to the line of sight 
direction 7. In the small angle approximation however, one assumes that k± is much larger compared to ku . We will denote 
the angle between the line of sight direction 7 and the wave vector k, by 6. The smoothing angle will be denoted by #o- Using 
the definitions that we have introduced above we can now express the smoothed projected two-point correlation function 
(Limber 1954, Peebles 1980, Kaiser 1992, Kaiser 1998): 

(«(7i)«(7 2 ))c = r^^M I 4k ex P(^) HzLWW)- (4) 
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Where we have introduced a new notation 1 = r(x)kx which denotes the scaled wave vector projected on the surface of the 
sky. The average of the two-point correlation function (k 2 ) smoothed over an angle do with a top-hat smoothing window 
W2(Wo) is useful to quantify the fluctuations in k s which is often used to reconstruct the matter power spectrum -P(k) (Jain 
& Selzak 1997). 
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Using the standard perturbative techniques, it is possible to compute the normalized cumulants or % parameters, for the 
convergence field in the large smoothing angle regime. 

2.1 The Probability Distribution Function of k(6o) 

To construct the probability distribution function (pdf) it is necessary to have the sn parameters for the convergence maps 
computed to an arbitrary order. These computations have already been done by Bernardeau et al. (1996) and Waerbeke et 
al. (1998). 

{k 2 (i)) = C2[I k J (6) 

< K 3 ( 7 )> = S rC3[^ o ] (7) 

<^ iv (7)) = .4 T C J v[/^ 1 ] (8) 

Where the sjv parameters [f] correspond to the case of 3D gravitational dynamics with 2D perturbations (Bernardeau 1995), 

T Throughout this paper we will use lower case letters for quantities related to both the convergence k and its reduced counterpart 
r\. We will show that the statistics associated with the reduced convergence map is very similar 
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and Ct denotes the line of sight integration. 
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The spectral index n is however same as that of the three dimensional (local) power spectral index. These results have already 
been tested using numerical simulation by Munshi et al. (1999). Typically for a degree scale smoothing, we will use n = —1.3. 
The projection effects are incorporated in the line of sight integration. Such a separation of dynamical part and the geometrical 
part occurs only when we replace the actual initial power spectra with the local power law spectra (Bernardeau, 1995). 
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Finally we can write down the normalized cumulants for the convergence field as: 
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The generating function approach developed by Balian & Schaeffer (1989) and later extended by Bernardeau & Schaeffer 
(1992) was used to compute the whole hierarchy of sm parameters for the three dimensional matter distribution and the 
projected matter distribution (Bernardeau 1995). The generating function is related to the vertices which appear in the tree 
representation of the correlation hierarchy which appear in the quasi-linear evolution of gravitational clustering. 
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Such a hierarchy also appears in the highly non-linear regime but with different tree amplitudes. The generating function for 
the sjv parameters can similarly be represented as (Balian & Schaeffer 1989): 
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The function <f> (y) satisfies the constraint si = S2 = 1 necessary for the normalization of the PDF. These two generating 
function are known to be related to each other by the following set of equations (Balian & Schaeffer 1989 , Bernardeau & 
Schaeffer 1992): 
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Hence once we know the generating function G(t) it is possible to compute the whole hierarchy of the projected sn parameters. 
It was shown that the generating function takes a particularly simple form (Bernardeau 1995): 
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with K a — (\/l3 — l)/2 for 2D perturbative dynamics of gravitational clustering in a 3D universe. This particular result hold 
for unsmoothed density field but it is possible to incorporate the effect of smoothing into account if we assume that the 
smoothing window is a top-hat function (Bernardeau 1994, Bernardeau 1995). 
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For the variance a at length scale 7?o will use the local power spectral index n so that a(R ) oc R 
The VPF and the PDF can be related to each other by the following equation (Balian & Schaeffer 1989): 
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to the projected two dimensional density fields. The Moments of these fields will be represented by a super-script PT to denote their 
perturbative origin, where as moments associated with the convergence map k will be denoted by the superscript n. 




Figure 1. The Probability distribution function associated with different realistic scenarios. The smoothing angle is fixed at 9q = 1°. 
The source is placed at z s = 1 for each of these models. The slope of correlation function 7 = 1.7 was assumed. Dots represent results 
from the log- normal fitting form with same variance and fc m in- 



The above expressions are clearly suitable for the density field and similar calculations can be done for convergence maps 
and the generating functions for the cumulants of the convergence maps can then directly be related to its counterpart for 
the density field. The details of such calculations can be found in Munshi & Jain (1999) and Valageas( 1999b). Our notations 
follow that of Munshi & Jain (1999) more closely. These studies clarified that it is useful to define a reduced convergence map 
r\ smoothed with an angle 9q; 

^=2^=1+^, (21) 
where the minimum value of the convergence along any line of sight direction can be defined by the following equation: 

Kmin = ~ I d,Xu(x)- ( 22 ) 
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Under certain approximation r\ will have the exactly same one-point statistics as the density field 1 + 5. 

$v(y) = <t> PT (y) (23) 

Notice that the quantity k min does not depend on the smoothing angle 9q and hence will only depend on the source red-shift 
z s and the definition of reduced convergence remains same both in the quasi-linear and the highly non-linear regime and is 
not connected with the convergence or divergence of the perturbative series. 

In general the expression connecting the cumulant generator for r\ ,i.e. & v can be related to 4> by the following approxi- 
mation. 
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Clearly this expression does not depend on whether we are in the quasi-linear regime or in the highly non-linear regime. 
Only changes which are to be made is that of <j>. In quasi-linear regime the quasi-linear phi can be computed from quasi-linear 
generating function for tree vertices defined earlier and in highly non-linear regime it needs to replaced with its non-linear 
counterpart. The simplification mentioned above can be achieved by replacing the red-shift dependence of different integrands 
with their values at a median red-shift. Numerical computations seems to suggest that such an approximation works remarkably 
well (Valageas 1999b, Munshi & Jain 1999a). 

It is possible to have different approximate asymptotic expressions for different limiting values of r\. However here we 
have performed a direct numerical integration and results for different cosmologies and various smoothing angles are presented 
in the Figure-1 and Figure-2. The different curves represent various versions of (standard, lambda, omega and r) cold dark 
matter models. The details of various parameters charecterising these dark matter models can be found in Munshi & Jain 
(1999a). 




Figure 2. The Probability distribution function associated with different realistic scenarios. The smoothing angle is fixed at do = .5° 
in the left panel and at 80 = .25° in the right panel . The source is placed at z B = 1 for each of these models. The slope of correlation 
function 7 = 1.7 was assumed. Dots represent results from the log- normal fitting formula with same variance and k m i n . 



As in the case of quasi- linear regime we now get: 



53 — S3 / K m i n (25) 

54 = S 4 / K min (26) 
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Although these expressions are very similar to their small angle counterparts (Munshi & Jain 1999a, Valageas 1999), however 
it is not very difficult to notice that there is a sharp difference between the two. For larger smoothing angle the statistics 
of reduced convergence map follows that of projected density field which in turn is related to the results associated with 
the two dimensional perturbations in a three dimensional universe. However in case of very small smoothing angles when 
we are probing the non-linear regime, the statistics of the convergence maps is that of full three dimensional density field. 
The intermediate regime is more difficult to model and it interpolates these two regimes in a smooth manner. Numerical 
investigations of S3 was done by Gaztanaga & Bernardeau (1998) for SCDM model for various smoothing angle. We find a 
good agreement between our results and their numerical measurements for angular scales as small as 8 = 15'. Which means 
that our analytical results for pdf and bias (based on tree-level calculations), may even be valid for much smaller angular 
scales. A detailed comparison however will require a large simulation patch and is not available at present. 



2.2 The bias associated with the convergence field k(6q) 

While it is clear that one-point studies of convergence field is interesting it was shown by Munshi (2001) that such results can 
also be obtained for the case of two-point quantities such as bias associated with the convergence map 6(«) and its moment, 
the cumulant correlators c pq (see Munshi, Melott & Coles 1999 for detailed descriptions of generalized cumulant correlators 
based on various hierarchical models). 

K<2)dK\d,K,2 = p(Kl)p(K2)(l + &( K l)£l2&( K 2))dKlGk2 (28) 

We have used shorthand notations k\ = 1(71), ^2 = ^(72) and £12 = (^(71)^(72)) to simplify the above expression. A similar 
expression can be obtained for the reduced convergence map r\. 

Our studies in the highly non-linear regime have shown that the hierarchical ansatz can be used to understand how the 
'hot-spots' in convergence maps are correlated and how such correlations can be used to study the similar correlation for the 
underlying matter distribution and hence the bias associated with the collapsed objects. We will show that a similar study 
is also possible in the quasi-linear regime with appropriate change in the generating function as was the case for non-linear 
regime. 



The normalized cumulant correlators or CpJ parameters have been calculated for dark matter density perturbations using 
direct perturbative calculations for low orders and for arbitrary orders using the generating function techniques. We can define 
these parameters by the following relation 

<«S(7iK(7 2 )>c == c^>? +9 - 2) <^( 7l )M7 2 )>. (29) 

These parameters were computed by Bernardeau(1996) for the density field and were found to obey a generic factorization rule 
c P q = Cpic q i for the case of 3D perturbations. Similar computations were later generalized for the case of 2D perturbations. 
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For projected fields we will get a prefactor resulting from the line of sight integration as was the case for the one-point 
cumulants. 
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where we have introduced a new quantity I K g 12 which incorporates the line of sight integration effects for the two-point 
quantities, such as the cumulant correlators and bias which we are considering here. 
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ft can be shown that the role of the generating function for the normalized cumulant correlators for the (unsmoothed) density 
field 1 + S Cpi is played by the quantity r itself. 
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In a very similar way we can define the generating function c pq for the convergence map n by the following relation: 
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However such a factorization is not possible for the projected density fields, without some simplifying assumptions and in 
general we have to use the generating function for CpJ which we denote as (3 K . Using relations we have already presented 
above we can finally write: 
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As was the case with one point cumulants we can use the reduced convergence ij(d ) to simplify the analysis because as 
pointed out earlier the statistics of reduced convergence map is very similar to that of the underlying projected density field. 
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If we use the same techniques which we have used to simplify the integrals associated with one the point cumulants (see 
Munshi & Jain 1999a, and Munshi 2000) where we replaced the red-shift dependence in integrands with their values at a 
median red-shift we can write: 



#,(2/1,2/2) = PMPM = /3 pt (j/i)/? pt M 



(38) 




Figure 3. Bias associated with different realistic scenarios. The smoothing angle is fixed at 80 = 1°. The source plane is placed at z s = 1 
for each of these models. The slope of the correlation function 7 = 1.7 was assumed. Dots represent predictions from log-normal model. 



This is exactly the same result which we have obtained for small smoothing angle, however results we have derived here are 
valid for large smoothing angles and hence the generating function (3 FT (y) corresponds to the case of perturbative regime. If 
we neglect the effect of smoothing the generating function f3 PT (y) is simply the function r(y) (Bernardeau 1996) however for 
smoothed convergence maps we can write (Bernardeau 1996): 
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Finally the bias b(rj) can be expressed as: 
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The generating function for r\ and projected density field are same under certain simplifying approximation, it shows that the 
bias associated with the reduced convergence is very similar to the statistics of the projected density field. This is in direct 
contrast with what was obtained using the nonlinear theory for small smoothing angle where the statistics of hot-spots in 
projected map are more closely related to the underlying three dimensional density field and not to the projected density 
field. Finally the bias associated with the convergence k can be expressed by the following equation. 
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Clearly it is possible to integrate the exact equation which represents joint moments exactly and hence compute the exact 
bias function for large smoothing angle from /3(y), but it was found by Munshi(2000) that simplification of replacing all the 
integrals by their approximate values, works remarkably well, when tested against numerical simulations. We will adopt the 
same technique for this study. Result of such numerical integration for different cosmologies are presented in Figure-3 and 
Figure-4. Instead of plotting the differential bias function we have plotted its integral part which simply represent the bias 
associated with spots which cross certain threshold in reduced convergence. Such a bias is easier to measure from numerical 
maps. We find that the bias predictions based on log-normal models also match perturbative calculations reasonably well. 

While in the non-linear regime it is possible to extend such studies to even higher order where one gets the complete picture 
for the many-point statistics of over-dense cells, similar analysis is not possible in quasi-linear regime as the tree-vertices v n 
defined above are dependent on shape factors. 

Different sources of noise affect the measurement of one-point quantities such as pdf and two-point quantities such as 
bias from observational data or numerical simulations. While for small angular scales the major concern is always the noise 
due to intrinsic ellipticities of galaxies at larger smoothing angle it is finite volume corrections which starts playing the major 
role. It manifests itself mainly on large n tail of the pdf and bias and depends on the presence or absence of a rare cluster in 
the catalog. The large k, tail generally shows large fluctuations before showing an abrupt cutoff for a certain threshold beyond 
which there are no more over-dense cells in the catalog. Similar effects are observed in the determination of bias too. 




Figure 4. The bias function associated with different realistic scenarios. The smoothing angle is fixed at 9q = .5° in the left panel and 
at 60 = .25° in the right panel . The source plane is placed at z s = 1 for each of these models. The slope of correlation function 7 = 1.7 
was assumed. Dots represent results from the log- normal fitting formula with the same variance and k m i n . 



3 DISCUSSION 

Perturbative calculations have already been used to compute the lower order hierarchy of sjv parameters associated with the 
convergence maps in quasi-linear regime. On the other hand hierarchical ansatz was shown to be useful in highly non- linear 
regime. In this particular paper we have shown how the recent development in highly non-linear regime can also be used 
to study the quasi-linear regime by introducing a reduced convergence field to construct the full pdf and bias associated 
with convergence maps in the quasi-linear regime. We have shown that a reduced convergence field in fact is related directly 
with the two-dimensional projected density field. It is known from earlier studies that the generating function corresponding 
to the two-dimensional density field is connected with the two dimensional spherical collapse. Which indicates that such a 
generating function now can be used to predict the pdf and bias of smooth convergence maps in the quasi-linear regime. A 
quick comparison with corresponding results in highly non-linear regime shows that while statistics of projected density for 
small smoothing is intimately related with the full three dimensional density field at larger smoothing angle it is related to 
the projected two dimensional density field. 

From observational point of view our results will have direct significance. While most recent observational studies cover 
rather small patches of the sky necessary for quasi-linear approximations to be valid, it is hoped that larger and larger patches 
of the sky will be covered in future, where one can apply perturbative techniques. At smaller angular scales it is well known 
that noise due to intrinsic ellipticity of galaxies will play a major role in determination of the cosmological parameters, however 
such noise will be less significant for studies at larger angular scales where our results will have direct relevance. Although 
signal variance will start declining with larger and larger smoothing angle and direct determination of weak lensing signal at 
degree scale might seem difficult at present, in the intermediate length scale #0 = 15' where the hierarchical ansatz starts to 
break down our analysis using perturbative techniques will provide much needed theoretical insight. 

Most of the presently available numerical maps for cosmic shear are generated from high resolution numerical simulation 
and are generally correspond to few degree patches in the sky which are suitable for testing analytical predictions for smaller 
smoothing angles, and recent comparison of various statistics were found to be very accurately reproduced by non-linear 
theory. Such studies for larger smoothing angles will have to wait till much larger numerical maps are available. There have 
been some attempts in this direction using higher order Lagrangian theories which replaces the exact gravitational dynamics. 
Such models although less accurate can reproduce the basic features of convergence maps in larger smoothing angles and 
are definitely very cost effective. For reproducing the statistics of convergence maps made using higher order Lagrangian 
perturbation theory we just need to change the exact generating function which we have used in this study by corresponding 
smoothed generating function for perturbative Lagrangian dynamics for the case of projected perturbations and it can then 
be used to compare with numerical studies. It can be shown that the frozen flow approximation as proposed by Matarrese 
et al(1992). will produce a log-normal pdf in the quasi-linear regime (Munshi, Sahni & Starobinsky 1994), hence such an 
approximation can also be used to generate convergence maps. 



Results presented here can be generalized very easily for the case of projected galaxy survey. A detailed study will be 
presented elsewhere. 
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APPENDIX: LOG-NORMAL DISTRIBUTION 



Lognormal distribution provides a phenomological description of one -point probability ditribution function (see Coles & Jones 
1991) and its generalisation to compute the bias function associated with overdense region (Taruya et al. 2001). Although 
inherently local in nature we found that it can also be used as a good fitting function for both one-point probability distribution 
function (see Bernardeau & Kofman 1995 for similar comparison in 3D for density PDF) and the bias function associated 
with the convergence maps. For the projected density contrast S (which is related to the reduced cnvergence field r\ = 1 + ft) 
we can write the pdf as: 



P(5)dS ■ 



V2ttE 2 



exp 



2 



2£ 



(1+5) 



d5. 



Similarly the joint probability distribution function can be written as (Kayo et al.2001): 



P{S 1 ,S 2 )dS 1 dS 2 



2ir y / (T, 2 - X 2 ) 



exp 



E(A? + Al) 



2AAiA 2 



2(£ 2 -X 2 ) J + (1 + &) 



-dSi dS 2 



Where we have introduced the following notations: 

X = \n{l + U), £=(l + ff„j), A, = In [(1 + 5i)V (1 + &ni)\ ■ 



(42) 



(43) 



(44) 



£ n i and (,u n denotes the non-linear and linear correlation functions and can be computed using suitable modelling of non-linear 
evolution e.g. Peacock & Dodds (1996). To compute the one-point and two-point pdf's for k from above expressions we have 
to apply the transformation rj — > n (which involve computation of K, m in) discussed before. A detailed comparison of bias from 
log-normal model and from hierarchical ansatz for the case of density field will be presented elsewhere. 



